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Abstract 

This work presents the application of the power series method (PSM) to find solutions of partial differential-algebraic 
equations (PDAEs). Two systems of index-one and index-three are solved to show that PSM can provide analytical 
solutions of PDAEs in convergent series form. What is more, we present the post-treatment of the power series 
solutions with the Laplace-Pade (LP) resummation method as a useful strategy to find exact solutions. The main 
advantage of the proposed methodology is that the procedure is based on a few straightforward steps and it does 
not generate secular terms or depends of a perturbation parameter. 
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Introduction 

As widely known, the importance of research on par- 
tial differential-algebraic equations (PDAEs) is that many 
phenomena, practical or theoretical, can be easily mod- 
elled by such equations. Those kinds of equations arise in 
fields like: nanoelectronics (Bartel and Pulch 2006), elec- 
trical networks (Ali et al. 2005, 2003; Gunther 2000) and 
mechanical systems (Simeon 1996), among others. 

In recent years, PDAEs have received much attention, 
nevertheless the theory in this field is still young. For 
linear PDAEs the convergence of Runge-Kutta method 
is investigated in (Strehmel and Debrabant 2005). The 
numerical solution of linear PDAEs with constant coeffi- 
cients and the study of indices are given in (Lucht et al. 
1997a, 1997b; Lucht and Strehmel 1998; Lucht et al. 1999). 
Linear and nonlinear PDAEs are characterized by means 
of indices which play an important role in the treatment 
of these equations. The differentiation index is defined 
as the minimum number of times that all or part of the 
PDAE must be differentiated with respect to time, in order 
to obtain the time derivative of the solution, as a con- 
tinuous function of the solution and its space derivatives 
(Martinson and Barton 2000). 
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Higher-index PDAEs (differentiation index greater than 
one) are known to be difficult to treat even numerically. 
Often such problems are first transformed to index-one 
systems before applying numerical integration methods. 
This procedure called index-reduction, can be very expen- 
sive and may change the properties of the solution. Since 
applications problems in science and engineering often 
lead to higher-index PDAEs, new techniques are required 
to solve these problems efficiently. 

Modern methods like homotopy perturbation method 
(HPM) (He 1999, 2000, 1998; Vazquez-Leal et al. 2012), 
homotopy analysis method (HAM) (Guerrero et al. 
2013), variational iteration method (VIM) (Khan Y, et al. 

2012) , generalized homotopy method (Vazquez-Leal 

2013) , among others, are powerful tools to approxi- 
mate nonlinear and linear problems. The HPM has been 
successfully applied to solve various kinds of nonlinear 
problems in science and engineering, including Volterra's 
integro-differential equation (El-Shahed 2005), nonlin- 
ear differential equations (He 1998), nonlinear oscilla- 
tors (He 2004), partial differential equations (PDEs) (He 
2005a), bifurcation of nonlinear problems (He 2005b) and 
boundary-value problems (He 2006). Recently, the mod- 
ifications of the HPM have been used to solve DAEs 
(Aminikhah and Hemmatnezhad 2011; Asadi et al. 2012; 
Salehi et al. 2012; Soltanian et al. 2010). Nevertheless, 
the power series method (PSM) (Forsyth 1906; Ince 1956) 
is a well-known classic straightforward procedure from 



Springer 



©2014 Benhammouda and Vazquez-Leal; licensee Springer. This is an Open Access article distributed under the terms of the 
Creative Commons Attribution License {http://creativecommons.0rg/licenses/by/2.O), which permits unrestricted use, 
distribution, and reproduction in any medium, provided the original work is properly credited. 



Benhammouda and Vazquez-Leal SpringerPlus 2014, 3:137 
http://www.springerplus.eom/content/3/1/137 



Page 2 of 9 



literature that can be applied successfully to solve differen- 
tial equations of different kind: linear ordinary differential 
equations (ODEs) (Coddington 1989; Forsyth 1906; Ince 
1956; Kreyszig 1999), nonlinear ODEs (Biazar et al. 2005; 
Fairen et al. 1988; Filipich and Rosales 2002; Filipich et al. 
2004; Guzel and Bayram 2005; Kreyszig 1999) and linear 
PDEs (Kurulay and Bayram 2009), among others. This 
method establishes that the solution of a differential 
equation can be expressed as a power series of the inde- 
pendent variable. 

In this paper we present the application of a hybrid 
technique combining PSM, Laplace Transform (LT) and 
Pade Approximant (PA) (Barker 1975) to find analytical 
solutions for PDAEs (Ebaid 2011; Gokdogan et al. 2012; 
Merdan et al. 2011; Momani and Erturk 2008; Momani 
et al. 2009; Sweilam and Khader 2009; Tsai and Chen 
2010; Yamamoto et al. 2002). Solutions to PDAEs are 
first obtained in convergent series form using the PSM. 
To improve the solution obtained from PSM's truncated 
series, we apply LT to it, then convert the transformed 
series into a meromorphic function by forming its PA. 
Finally, we take the inverse LT of the PA to obtain 
the analytical solution. This hybrid method (LPPSM), 
which combines PSM with Laplace-Pade post-treatment 
greatly improves PSM's truncated series solutions in 
convergence rate. In fact, the Laplace-Pade resumma- 
tion method enlarges the domain of convergence of the 
truncated power series and often leads to the exact 
solution. 

It is important to remark that LPPSM can obtain 
exact solutions without requiring the index-reduction 
of the PDAEs. The proposed method does not pro- 
duce noise terms also known as secular terms as 
the homotopy perturbation based techniques (Soltanian 
et al. 2010). This greatly reduces the volume of com- 
putation and improves the efficiency of the method in 
comparison to the perturbation based methods. What is 
more, LPPSM does not require a perturbation param- 
eter as the perturbation based techniques including 
HPM. Finally, LPPSM is straightforward and can be pro- 
grammed using computer algebra packages like Maple or 
Mathematica. 

The rest of this paper is organized as follows. In the 
next section we illustrate the basic concept of the PSM. 
The main idea behind the Pade approximant is given 
in section "Pade approximant". In section "Laplace-Pade 
resummation method", we give the basic concept of 
the Laplace-Pade resummation method. The application 
of PSM to solve PDAE systems is depicted in section 
"Application of PSM to solve PDAE systems". In 
section "Test problems", we apply LPPSM to solve 
two PDAEs problems of index-one and index-three. In 
section "Discussion", we give a brief discussion. Finally, a 
conclusion is drawn in the last section. 



Basic concept of power series method 

It can be considered that a nonlinear differential equation 
can be expressed as 



A(u) -f{t) =0, teQ,, 
having as boundary condition 



ih) ) = 



t € r, 



(1) 



(2) 



where A is a general differential operator, f(t) is a known 
analytic function, B is a boundary operator, and T is the 
boundary of domain Q. 

PSM (Forsyth 1906; Ince 1956) establishes that the solu- 
tion of a differential equation can be written as 



u (t) = ^ u n t", 

n=0 



(3) 



where Wn, U\, . . . are unknowns to be determined by series 
method. 

The basic process of series method can be described as: 

1. Equation (3) is substituted into (1), then we regroup 
the equation in terms of powers of t, 

2. We equate each coefficient of the resulting 
polynomial to zero. 

3. The boundary conditions of (1) are substituted into 
(3) to generate an algebraic equation for each 
boundary condition. 

4. Aforementioned steps generate an algebraic linear 
system for the unknowns of (3). 

5. Finally, we solve the algebraic linear system to obtain 
the coefficients uq, U\, . . . 

Pade approximant 

Given an analytical function u{t) with Maclaurin's expan- 
sion 



u(t) = ^2 u„t", 0 < t < T. 



(4) 



n=0 



The Pade approximant to u (£) of order [L, M] which we 
denote by [L/M] u (t) is defined by (Barker 1975) 



[L/M\ u (t) 



Po+Pit+ ...+p L t L 
l + qit+ ... + q M t M ' 



(5) 



where we considered qo = 1, and the numerator and 
denominator have no common factors. 

The numerator and the denominator in (5) are con- 
structed so that u (t) and [L/M] u (t) and their derivatives 
agree at t = 0 up to L + M. That is 



u(t)-[L/M] u (t) = 0{t L+M+1 ). 



(6) 
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From (6), we have 

M L 

U (t) J2 Int" ~ Y^Pnf = O {t L+M+l ) . (7) 

(1=0 M=0 

From (7), we get the following algebraic linear systems 

ulIi + ■ ■ ■ + Ui-M+llM = —ul+i 
UL+iqi + ■ ■ ■ + ul-m+iVm = —ul+2 
■ . (8) 

. UL+M-lH + ■ ■ • + UiqM = -ul+m, 

and 

p 0 = u 0 

pi = u\ + u 0 qi 

. (9) 
Pl = u L + u L -\q\ + . . . + u 0 qL- 

From (8), we calculate first all the coefficients q„, 1 < n < 
M. Then, we determine the coefficients p„, 0 < n < L 
from (9). 

Note that for a fixed value of L + M + 1, the error (6) is 
smallest when the numerator and denominator of (5) have 
the same degree or when the numerator has degree one 
higher than the denominator. 

Laplace-Pade resummation method 

Several approximate methods provide power series solu- 
tions (polynomial). Nevertheless, sometimes, this type of 
solutions lacks of large domains of convergence. There- 
fore, Laplace-Pade (Ebaid 2011; Gokdogan et al. 2012; 
Merdan et al. 2011; Momani and Erturk 2008; Momani 
et al. 2009; Sweilam and Khader 2009; Tsai and Chen 2010; 
Yamamoto et al. 2002) resummation method is used in lit- 
erature to enlarge the domain of convergence of solutions 
or inclusive to find exact solutions. 
The Laplace-Pade method can be explained as follows: 

1. First, Laplace transformation is applied to power 
series (3). 

2. Next, s is substituted by 1 /t in the resulting equation. 

3. After that, we convert the transformed series into a 
meromorphic function by forming its Pade 
approximant of order [L/M].L and M are arbitrarily 
chosen, but they should be of smaller value than the 
order of the power series. In this step, the Pade 
approximant extends the domain of the truncated 
series solution to obtain better accuracy and 
convergence. 

4. Then, t is substituted by 1 /s. 

5. Finally, by using the inverse Laplace 5 transformation, 
we obtain the exact or approximate solution. 



Application of PSM to solve PDAE systems 

Since many application problems in science and engi- 
neering are often modelled by semi-explicit PDAEs, we 
consider therefore the following class of PDAEs 

wit = <P (u, u x , u xx ) , (10) 
0 = f (u,u x ,u xx ),(t,x) e (0,T) x (a,b), (11) 

where u k : [0, T] x [ a, b] -+ R m <s k = 1, 2 and b > a. 
System (10)-(11) is subject to the initial condition 

u\ (Q,x) =g(x),a < x < b, (12) 

and some suitable boundary conditions 

B (u (t, a) , u (t, b) , u x (t, a) , u x (t, b)) = 0, 0 < t < T, 

(13) 

where g (x) is a given function. 

We assume that the solution to initial boundary- 
value problem (10)-(13) exists, is unique and sufficiently 
smooth. 

To simplify the exposition of the PSM, we integrate first 
equation (10) with respect to t and use the initial condition 
(12) to obtain 

f 

u\(t,x) - g(x) - j 4> (w, u x , u xx ) dt = 0. (14) 
o 

It is important to note that the time integration of 
equation (10) is not relevant to the solution procedure 
presented here, so one can apply the PSM directly to (10). 

In view of PSM, we assume the solution components 
Uk(t,x) ,k = 1, 2 to have the form 

u k (t, x) = a k>0 (x) + a kA (x) t + a ki2 (x)t 2 + (15) 

where a ki „ (x), k = 1, 2; n = 0, 1, 2, . . . are unknown 
functions to be determined later on by the PSM. 

Then substitute (15) into system (11)-(14) and equate 
the coefficients of powers of t in the resulting polyno- 
mial equations to zero to get an algebraic linear system for 
these coefficients. Finally, we use equation (15) to obtain 
the exact solution components u k , k = 1,2 as series. 
The solutions series obtained from PSM may have limited 
regions of convergence, even if we take a large number 
of terms. Therefore, we apply the Laplace-Pade resum- 
mation method to PSM truncated series to enlarge the 
convergence region as depicted in the next section. 

Test problems 

In this section, we will demonstrate the effectiveness and 
accuracy of the LPPSM presented in the previous section 
through two PDAE systems of index-one and index-three. 
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Nonlinear index-one system: Then, we substitute (22) into equations (17) and (21) to 

Consider the following nonlinear index-one PDAE which get 

arises as a similarity reduction of Navier-Stokes equations t 

(Budd et al. 1994) g ^ {%) f _ C0S7TX _ ff- ^ (x) t « dt 

M=0 J «- n 



">- 2 f^ (I6) + ((±^A(£^,A 

J Q \„=0 I \„=0 / 

- f feairf. (*)«") * 

{ \n=0 J 



U\t = Ui xx - U2U\ X + 

0 = Uix - U\, (17) 



dt 



where 0 < x < 1 and t > 0. 

System (16)-(17) is subject to the following initial condi 
tion 



o 

+ 2 // i^^hnix)^ dxdt = 0, (23) 



U\ {0, X) = COS TTX, 0<X<1, (18) 0 0 

oo 

and boundary conditions { a 2,n (*) ~ a i,n(x)) t" = 0, (24) 

H=0 

, „ s , „ s , , „ s „ where (') denotes the ordinary derivative with respect to 

u lx (t,0) = u lx (t,l) = u 2 (t,0) = u 2 (t,l) = 0,t>0. w 7 F 

(19) 

Equating the coefficients of powers of t to zero in (24) 
then solving the resulting equation for (X2,n (*) and using 
The exact solution of problem (16)-(19) is the boundary conditions (19), we have 

X 

u x (t,x) = e-* 2 * cos it x, a 2 ,„(x)= fa 1>n (x)dx, n = 0,1,2,... (25) 

u 2 (t,x) = (l/n)e~ n2t smitx, 0 < x < 1, t > 0. 

Now equation (23) can be written as a series 

,.rr . . r ■ I -, r,s , («1.0 (*) _ COS 7tx) + > I (Xl, n (x) — (l/«) 0^1 „ i(x) 

Since one time differentiation of equation (17) deter- v ' 



n-l 



mines u 2t in terms of u and its space derivatives, then 

PDAE (16)-(17) is index-one. Note that no initial condi- _ (i/«)^^ H (x) J t n = 0, 

tion is prescribed for the variable u 2 as this is determined k=0 / 

by the PDAE. (26) 

In order to simplify the exposition of the PSM presented 
in section "Application of PSM to solve PDAE systems" where 

to solve (16)-(17), we first integrate equation (16) with r 
respect to t and use the initial condition (18) to get & k -» (x) = ai * (x) (*) - a'i,„-i-* (*) J «U (*) dx 

o 

/ L xx - u 2 u lx + u\-2 ju\dx\ dt = 0. ~ 2 / aU iX) ai - n - 1 ~ k (X) ^ 

0 V 0 / 

(2i) Equating all coefficients of powers of t to zero in (26), 
yields ai,o 0*0 = cos irx and the recursive formula for 
a\,„ (x) 

n-l 

«i,„ (*) = (l/«) «!,„_! (*) + (l/«) ^ (*) , 



Ml (t,X) — COS JTX — / «i„ — «2«lx + 



In view of the PSM, we assume the solution components 
u k ,k = 1, 2 to have the form 



u k (t, x) = a k?0 (x) + a kt i (x) t + a ki2 (x)t 2 + (22) 



( 27 ) 



n = 1,2,3, ... 



where a^„ (x),k = 1,2; n = 0, 1,2, ... are unknown From recursion (27), we get ai,i (x) = -it 2 cos nx and 
functions to be determined later on by the PSM. «i,2 (*) = (^ 4 / 2 ) cos nx. 
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From equation (25), we get 0/2,0 (x) = (I /ti) sin tix, 
a 2,i (x) = —ti sin Ttx and 0:2,2 (*) = (^ 3 /2) sin 71 x Using 
(22) and the coefficients recently obtained, we have 



Uy (t,X) = ^1 - JT 2 t + i (-7T 2 i) 2 ^ 



COS 7TX 



(28) 



and 



Ul 



(t,x) = (l - JT 2 i+ i (-tt 2 ?) 2 ^ (l/7T)sin7T^. (29) 

Similarly, the coefficients ai )W (#) and a% n (x) for « > 3 
can be found from (27) and (25) respectively. 

The solutions series obtained from the PSM may have 
limited regions of convergence, even if we take a large 
number of terms. Accuracy can be increased by applying 
the Laplace-Pade post-treatment. First, we apply i-Laplace 
transform to (28) and (29). Then, we substitute s by 1/t 
and apply i-Pade approximant to the transformed series. 
Finally, we substitute t by 1/s and apply the inverse Laplace 
5-transform to the resulting expressions to get the approx- 
imate or exact solutions. 

Applying Laplace transforms to u\{t,x) and Ui(t,x) 
yields 



£[«i (f,*)] = (--— + — ) 



COS TIX, 



(30) 



and 



,1 Tl 2 Tl 4 \ 

L [u 2 (t,x)] = ( T + — ) (1/7T) sin ttx (31) 



For the sake of simplicity let s= 1/t, then 

C [ui (t,x)] = (t - Ti 2 t 2 + Ti 4 t 3 ) cos ttx (32) 

and 

C [u 2 (t,x)] = (t - Ti 2 t 2 + tt 4 ? 3 ) (1/tt) sin ttx (33) 

All of the [L/M] i-Pade approximants of (32) and (33) with 
L > 1 and M > 1 and L + M < 3 yield 



[L/M] U1 (t,x) 



1+Tl 2 t 



COS TlX, 



and 



[L/M] 



ft*)=( r ^)(l/ ff )A 



(34) 



(35) 



Now since t = 1/s, we obtain [L/M] U1 and [L/M] U2 in 
terms of s as follows 

[L/M] U1 {t,x) = {Ti 2 +s)~ cos tix, (36) 

[£/M]„ 2 (f,*) = (ti 2 +s)~ 1 (l/7T)sin7rx (37) 

Finally, applying the inverse LT to the Pade approximants 
(36) and (37), we obtain the approximate solution which is 
in this case the exact solution (20) in closed form. 



Linear index-three system: 

Consider the following index-three PDAE system 

= u\ xx + W3 sin tix, (38) 

«2tt = U2xx + «3 COS TIX, (39) 

0 = U\ sin tix + U2 cos tix — e _t , (40) 

where 0 < x < 1 and t > 0. 

System (38)-(40) is subject to the following initial condi- 
tions 

u\(0,x) = smnx u\ t (Q,x) = — sinnx, (41) 

U2 (0,X) = COSTIX U2t (0, x) = — cos tix 0<x<l, 

(42) 

and the boundary conditions 

Mi (t,0) = Ml (t,l) = 0, 
M2 (t, 0) = -w 2 (t, 1) = e _t £ > 0. 

The exact solution of problem (38)-(43) is 

Mi (t, x) = e~ l sin tix U2 (t, x) = e _t cos ttx, 
m 3 (t,x) = (l + ?r 2 ) e _t , 0 < x < 1, £ > 0. 



(43) 



(44) 



Since three time differentiations of equation (40) deter- 
mine k 3( in terms of the solution u and its space deriva- 
tives, then PDAE (38)-(40) is index-three. Therefore, this 
PDAE is difficult to solve numerically. Moreover no ini- 
tial condition is prescribed for the variable M3 as this is 
determined by the PDAE. 

In order to simplify the exposition of the LPPSM pre- 
sented in section "Application of PSM to solve PDAE sys- 
tems" to solve (38)-(43), we first integrate equations (38) 
and (39) twice with respect to t and use the initial condi- 
tions (41)-(42) to get 

t t 

Mi {t, x) — sin tix -\- 1 sin tix — j j u\ xx 

0 0 

+ M3 sin Tixdtdt = 0, (45) 
t t 

U2 (t, X) — COS TlX + t COS TIX — j j U2x X 

0 0 

+ M3 cos Tixdtdt = 0. (46) 



In view of the PSM, we assume the solution components 
Uk(t,x) ,k = 1,2,3 to have the form 

u k (t, x) = a kfi (x) + a kA (x) t + a k>2 (x)t 2 + (47) 

where a k>n (x) , k = 1, 2, 3; n = 0, 1, 2, . . . are unknown 
functions to be determined later on by the PSM. 
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«=0 



Substituting (47) into equations (40), (45) and (46) we and the nonsingular algebraic linear system for the 

get the system unknown functions ai ?n , oi2,n and a^, n -2 

oo * * oo , , a 3 ,„_2(*)sinjr* a" 2 (*) 

^ a\ >n (x) t" — sin nx + t sin jtx — / / y (x) fdtdt '" (n — Y)n (n-l)n' 

0 0"=° gjg-^W costt* _ <* 2 ',«-2 (*) 

a2 '" (X) (n-l)n ~ (n-l)n' (52) 

/ / Ta 3 ,„(x)t n dtdt = 0, (-I)" 

J J *—L a\ n (x) sin 7TX + «2 n (X) COS J!* = 

o o n =° «! 

(48) for n = 2, 3, . . . 

Solving system (52) exactly, we obtain the recursions 

oo * j. oo (-1)" . S„(x)COSnX 

^2 u%n (x) t" — cos nx + t cos 7M — / / a 2i „ Qs) t n dtdt a i,n(X) — ^ smnx+ (y l — \^ n 

„_n •/ •/ , — A - . „ . . 



t f 

sin jrx 



«=o 



* t 

COS 7TX 



(—1)" 5„ (x) sin J!* 

a 2;M (#) = cos 7i x - 



, . n\ (« - 1) « (53) 

/ / J2a 3 , n (x) fdtdt = 0, (-1)" „ 

o o H - u (« — 2) ! 

(49) — a 2 „_ 2 (#) cos 7r#, 

an d where 5„ (x) = c/'{ n _ 2 (x) cos 7T# — « 2 «-2 (*) s ^ n 7rx 

go For « = 2, 3, 4, we have S„ (x) = 0 and hence 

(ai,„ (#) sin + a 2 ,H (#) cos nx) t" — e~ t = 0, 2 1 

h=o a i,2 (x) = - sin ttcv, 0:2,2 (*) = - cos nx, 

(5°) ^1.2 
where (') denotes the ordinary derivative with respect to 

1 

ai 3 (x) = — sin n 
6 

«3,1 (*) = - (l + 7T 2 ) , 



X. 



1 1 

/™% , . . «13(*) = SHUT*, C(23(X) = COS7T#, 

System (48)-(50) can be rewritten as series 6 6 



(ai,o (*) — sin nx) + (ai,i (*) + sin jtx) t an( j 



00 / 



V- I <»- 2 W + «3,»-2 (x)sinTTX \ I 1 

"2^^ (h-1) h <xi,n(x)jt = 0, «i, 4 (3c) = — smjtx, a 2A {x)= — 



COS 



(«2,0 (x) — COS 7T;«) + (o!2,l (x) + COS TTje) t 



a 3 ,2 (*) = ^ (1 + JT 2 ) . 
Using (47) and the coefficients recently obtained, we get 



oo / „,// 



| «2,«-2 (X) + <*3,«-2 (X) COS ( x )\ t » = 0 / 111 

" («-!)« a 2,«WI , Ml (t,«) = (l-£ + -£ 2 Z: 3 + -Z: 4 I sin^, (54) 



K=2 \ 

(51) 



3! 4! 



( r\- , ^ a "2 (t, JC) = I 1 - t + ~t 2 - — t 3 + — t 4 | COS7T^, (55) 

yai : „ (x) sin nx + a 2 ,« (*) cos nx — J t" = 0. * v ' > ^ 1 2 3! 4! / 



(56) 



«=o 

and 

Equating the coefficient of powers of t to zero in (51) 
then solving the resulting system we find the coefficients M3 (f t x) = (l + n 2 ) ( 1 — t + -t 2 ) ■ 

a ki „(x), fork= l,2,3and« = 0,1,2, ... V 2 / 

Similarly, the coefficients a\ t „ (x), a 2 , n (x) and a3,„_2 (x) 
B y (*) = sinTr^, (*) = - sinTr^, for « > 5 can be found from (53). The solutions series 

obtained from the PSM may have limited regions of 
«2,0 (*) = costt^;, a2,i (*) = — cos7T^, convergence, even if we take a large number of terms. 
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Therefore, we apply the t-Pade approximation technique 
to these series to increase the convergence region. First 
t -Laplace transform is applied to (54), (55) and (56). 
Then, s is substituted by 1 ft and the t-Vade approximant is 
applied to the transformed series. Finally, t is substituted 
by 1 /s and the inverse Laplace s -transform is applied to 
the resulting expressions to get the approximate or exact 
solutions. 

Applying Laplace transforms to u\(t,x), U2(t,x) and 
U3 (t, x) yields 



C [ui (t,x)] = 

C [U2 (t,X)] = 

and 

£ [w 3 (t,x)] = 




For the sake of simplicity let 5 = I ft, then 

C [Ml (t,X)] = (t - 7t 2 t 2 + TT 4 ? 3 ) SHITT*, (60) 

C[u 2 {t,x)} = (t - 7t 2 t 2 + ?T 4 i 3 ) COS 71 X, (61) 

and 

C [w 3 (t,x)] = (1 + Jt 2 ) (t - it 2 t 2 + 7r 4 t 3 ) . (62) 

All of the [L/M] t-Pade approximants of (60), (61) and (62) 
with L > 1 and M>landl+M<3 yield 

[L/M] U1 (t,x) = (yY^j sin**, (63) 

[L/M] U2 (t,x) = (yY^j cositx, (64) 

and 

[L/M] U3 (t,x) = {l + it 2 )^^-^j. (65) 

Now since t = 1/s, we obtain [L/M] Ul , [L/M] U2 and 
[L/M] U3 in terms of 5 as follows 

[L/M] U1 {t,x) = {jT 2 +s)~ sin jtx, (66) 

[L/M] U2 (t,x) = (jt 2 costt*, (67) 

and 

[L/M] U3 (t,x) = (1 + it 2 ) (iz 2 + s)" 1 . (68) 

Finally, applying the inverse Laplace transform to the 
Pade approximants (66), (67) and (68), we obtain the 



approximate solution which is in this case the exact solu- 
tion (44) in closed form. 

Discussion 

In this paper we presented the power series method (PSM) 
as a useful analytical tool to solve partial differential- 
algebraic equations (PDAEs). Two PDAE problems of 
index-one and index-three were solved by this method 
leading to the exact solutions. The method has success- 
fully handled the index-three PDAE without the need for a 
preprocessing step of index-reduction. For each of the two 
problems solved here, the PSM transformed the PDAE 
into an easily solvable linear algebraic system for the coef- 
ficient functions of the power series solution. To improve 
the PSM solution, a Laplace-Pade (LP) post-treatement 
is applied to the PSM's truncated series leading to the 
exact solution. Additionally, the solution procedure does 
not involve unnecessary computation like that related to 
noise terms (Soltanian et al. 2010). This greatly reduces 
the volume of computation and improves the efficiency of 
the method. It should be noticed that the high complex- 
ity of these problems was effectively handled by LPPSM 
method due to the malleability of PSM and resumma- 
tion capability of Laplace-Pade. What is more, there is 
not any standard analytical or numerical methods to solve 
higher-index PDAEs, converting the LPPSM method into 
an attractive tool to solve such problems. 

On one hand, semi-analytic methods like HPM, HAM, 
VIM among others, require an initial approximation for 
the sought solutions and the computation of one or sev- 
eral adjustment parameters. If the initial approximation 
is properly chosen the results can be highly accurate, 
nonetheless, no general methods are available to choose 
such initial approximation. This issue motivates the use 
of adjustment parameters obtained by minimizing the 
least-squares error with respect to the numerical solution. 

On the other hand, PSM or LPPSM methods do not 
require any trial equation as requisite for the starting the 
method. What is more, PSM obtains its coefficients using 
an easy computable straightforward procedure that can 
be implemented into programs like Maple or Mathemat- 
ica. Finally, if the solution of the PDAE is not expressible 
in terms of known functions then the LP post-treatement 
will provide a larger domain of convergence. 

Conclusion 

This work presented LPPSM method as a combination 
of the classic PSM and a resummation method based 
on the Laplace transforms and Pade approximant. Firstly, 
the solutions of PDAEs are obtained in convergent series 
forms using PSM. Next, in order to enlarge the domain 
of convergence of the truncated power series, a post- 
treatment combining Laplace transform and Pade approx- 
imant is applied. This technique that we call LPPSM 
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greatly improves PSM's truncated series solutions in con- 
vergence rate, and often leads to the exact solution. Addi- 
tionally, PSM is an attractive tool, because it does not 
require of a perturbation parameter to work and it does 
not generate secular terms (noise terms) as other semi- 
analytical methods like HPM, HAM or VIM. 

By solving two problems, we presented the LPPSM as 
a handy tool with high potential to solve linear/nonlinear 
higher-index PDAEs. Additionally, the LPPSM does not 
require an index-reduction to solve higher-index PDAEs. 
Furthermore, we obtained successfully the exact solu- 
tions of such two problems highlighting the efficiency of 
LPPSM. What is more, the proposed method is based 
on a straightforward procedure, suitable for engineers. 
Finally, further research should be performed to solve 
other higher-index nonlinear PDAE systems. 
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